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Abstract 

<N . . 

J> I The crossing properties of the matrix elements of non-local operators, parameterized 

' by Generalized Parton Distribution, are considered. They are especially simple in terms 

. of the Double Distributions which are common for the various kinematical regions. As 

I a result. Double Distributions may be in principle extracted from the combined data 

in these regions by making use of the inverse Radon transform, known as a standard 
method in tomography. The ambiguities analogous to the ones for the vetor potential 
in the two-dimensional magneto-statics are outlined. The possible generalizations are 



d.cu.sed. 

^ I Introduction 

Deeply Virtual Compton Scattering (DVCS) [||, 0] is the cleanest hard process which is 
^ I sensitive to the Generalized Parton Distributions (GPD), the most popular version of which is 
' probably represented by the Skewed Parton Distributions (SPD), and has been the subject of 
extensive theoretical investigations for a few years. First experimental data became recently 
available (see e.g. ^) and much more data are expected from JLAB, DESY, and CERN 
in the near future. In the present note we shall limit ourself to the case of the pion target, 
although our consideration can be apphcable to any target. 

The crossing version is then provided by the process 7*7 tttt with a highly virtual 
photon but small hadronic invariant mass W. It was recently investigated in the framework 
of QCD factorization 0]- It allows to study the pion pair produced in the isoscalar 
channel, where the huge p-meson peak is absent. This process is analogous to the single 
pion production, described by the pion transition form-factor, being the long time object of 
QCD studies ^. In particular, the generalized distribution amplitude (GDA), describing 
the non-perturbative stage of this process, is a natural counterpart of the pion light cone 
distribution amplitude. From the other side, it may be considered as a crossing counterpart 
of pion SPD 1^ and related to it |TD| by making use of the suitable polynomial basis. 

The alternative (and quite elegant) description of non-perturbative stage of DVCS and 
other hard exclusive processes is provided by Double Distributions (DD) Jl], They 
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naturally explain the polynomiality of SPD and can be a good starting point for the model- 
building. At the same time, it is not clear how one can express DD in terms of the hadronic 
matrix element and extract them, at least in principle, from the data. 

In this note we study the crossing properties of the double distributions. We conclude, 
that both mentioned channels may be described by the common DD, while SPD and GDA 
may be obtained by its integration over the straight lines constituting the various angles 
with the coordinate axes. This would allow us to recover the DD from the SPD and GDA 
by making use of the inverse Radon transform - the mathematical tool widely used for 
computer tomography. We discuss the ambiguities resulting from the presence of the so- 
called Polyakov- Weiss (PW) terms in DD and find their close analogy with that for the 
recovering of the vector potential from the magnetic field. We also outline the possible 
apphcations of the method and its further development. 

Double distributions in the crossing related channels 

Let us start with the following symmetric representation of the matrix elements of twist-2 
non-local operators (c.f. |TU|) 

r-l-|x| 



J-l J\x\-l 



-1 J\x\-l 
r-1 rl~\x\ 



(p',-p|^(-|)7-#(|ri)|0) = {2P.z) j_^dx^^^^ dye-'^''^-'y''^/'f{x,y,A') 

+ {A-z) /' dx /' dye-'''^'-'y'^'/^g{x,y,A^y, 

J — 1 J\x\ — 1 



where we (in order to express the crossing properties in the most simple way) adopted the 
common notations for both channels: P = {p + p')/2,A = p' — p = 2^P + A^, while 
p,p' are the initial and final momenta of pion in the DVCS channel, respectively. In the 
7*7 — i> vrvf channel, A corresponds to the total momentum of pion pair while 2P is the 
relative momentum of pions. 

The skewedness parameter varies in the region < ^ < 1 for the DVCS channel and 
1^1 > 1 for the 7*7 — > tttt channel, respectively. In what follows we restrict ourselves to the 
leading twist level, neglecting the transverse momentum A_l. 

The crossing symmetry states that in the different channels the matrix elements are 
described by the common analytical functions, so that one may relate them by analytical 
continuation. One should note the crucial advantage of the DD for the studies of crossing 
properties, namely, that ^ does not enter explicitly to their definition. The only variable, 
which should be continued for DD is therefore A^ which is equal to t in the DVCS channel and 
to the squared invariant mass of the pion pair W"^ in the 7*7 tttt channel. The dependence 
on these variables is entirely non-perturbative and is a subject of separate investigation. This 
dependence is beyond the scope of our studies, and we shall put both t and W"^ to zero. This 
point is the common unphysical point for the both channels, so one may conclude, that 
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F{x, y, 0)) = /(x, y, 0); G{x, y, 0) = g{x, y, 0) 



(2) 



The Skewed Parton Distribution are related to the Double ones by the integration over 
straight lines [|TT|. The symmetric skewed distribution is related to symmetric double distri- 



butions F and G by the formula, being the straightforward generalization of |11 



H{z, = dx dy{F{x, y) + iG{x, y))5{z -x-iy). (3) 

In turn, the Generalized Distribution Amplitudes, describing the non-perturbative stage 
of hadron pair production, may be related to the respective Double Distributions in an 
analogous way: 

1/0 = dx Z'"'"^' dy{g{x, y) + - x/i - y) (4) 

Note the difference with the original definition of GDA P], where skewedness was defined 
by the fraction of the pair momentum carried by the one of the produced pions, while here 
it is defined (in the more symmetric way) by the difference of their momenta. As a result, 
the analytic continuation between $ and H is provided by the formula: 



i7(z,0 = 5^^n(0$(|,^) (5) 

Contrary to (^, SPD and GDA, which are depending on ^ explicitly, correspond to the 
different regions of this variable: ^ < 1 for SPD and \^\ > 1 for GDA. In what follows, we 
shall consider SPD in the 'extended' region, on the whole real axis, assuming that for |0 > 1 
it should be substituted by GDA from (|^), while the region — 1 < ,^ < is accessible by the 
use of symmetries. 

These discrete symmetries of the SPD and GDA may be now described in the unified 
way. Namely, the symmetries of DD with respect to its second argument 

F{x, y) = F{x, -y) ; G{x, y) = -G{x, -y), (6) 

correspond to either T-invariance for SPD, or charge conjugation invariance for GDA. This 
is not surprising, as the T-invariance, first studied in the case of the forward twist-3 matrix 
elements |jl2|, corresponds to the interchange of both hadrons and partons, in complete 
analogy with charge conjugation. Note that the reality of the scalar function (SPD or DD) 
plays the crucial role in such a derivation. At the same time, GDA may acquire the non- 
trivial imaginary phases due to cut in W^. However, it is only the real part of GDA which 
is described by crossing (^, while the imaginary parts requires the separate study of the 
analytical continuation in W^. Also, as soon as the chosen definition of GPD preserves P - 
invariance, such a similarity of T- and C- invariance is related to CPT-theorem. Note that the 
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action of the T transforms GDA to another objects: {0\ip y— ^ j 'j ■ zip y^n j \ —p',p), describing 

the appearance of hh pair in the initial, rather than in the final state. Unfortunately, they, 
strictly speaking, cannot be accessed in the process of /i/i-annihilation to the 7*7 pair, being 
the crossing version of DVCS and hh production, as the mass of virtual photon, which should 
provide the 'hard' scale, is constrained by the kinematics to be smaller than s = W^, which 
is assumed to be the 'soft' parameter. At the same time, one may consider the hh (in reality, 
probably pp) annihilation to the pair of two real photons at (moderately) large s as a crossing 
version of the large angle (real) Compton scattering, which is known to be described by GPD 
[^. This may give access to the pp GDA, very difficult otherwise. 

The relations (^ lead to the following symmetry properties for SPD and GDA 

H{x, = Hix, -0; Hz, = -H-z, -0, (7) 

The additional symmetry properties emerge if the pions are in the definite isospin state 
(where we assume that more general symmetry (^,0) also holds) 

F'='{x,y) = -F'=\-x,yy,G'=\x,y) = G^=°(-x,y), 
F'=\x,y) = F'=\-x,y);G'=\x,y) = -G'=\-x,y) 

H'='ix,0 = -H'='hx,0;^'=%z,0 = -$'=°(-^,e) = -<^'='{z,-0 

H'=\x,0 = H'=\-x,0;^'=\z,0 = ^'=\-z,0 = -$'='(^,-0- (8) 

For simplicity, we shall limit ourselves for the time being to the case of DD F{x, y) only, 
neglecting the Polyakov- Weiss [|TU] term, resulting from the function G{x,y). This approxi- 



mation is self-consistent in the case of / = 1 



Double distributions from the inverse Radon transform 

We now make the key observation, that the relation between SPD and DD is nothing else 



than a particular case of the Radon transform [0, |T5|, |T6|, 0. As soon as all the possible 
lines crossing the compact region are considered, one is dealing with the mapping between 
functions of two variables (because each line is characterized just by two real numbers). As 
soon as the transformed function is continuous, it may be recovered by making use of the 
inverse Radon transform. This procedure is known to be the key ingredient of the numerous 
applications [l^, covering, say, medicine, optics and geophysics. We are going to suggest its 
application in (non-perturbative) QCD. 

As soon as Radon transform is a new instrument in this field, it makes sense to present 
briefiy the derivation of its inversion, which is actually very simple. It is more convenient to 
present it by parameterizing the straight line in a slightly different way: by using the unit 
vector ^ = (cos0, sm0) orthogonal to it, instead of skewedness ^, and its distance from the 
origin p instead of the argument of the skewed distribution z. The correspondence to (^ is 
obvious. 
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R{p, = f dx dyfix, y)S{p - ff), (9) 

^ = tg<P,z = p/cos<j);H{z,0 = R{p,0\cos<P\. (10) 



The most convenient way to invert the transform is to use its relation to Fourier transform. 
Indeed, the Fourier transform of the function f{x) = f{x,y) may be written in the form 



F{q)= yrfW"V(^) = J dtS{t-xq) J xe''"' f (x) (11) 

Here the integration is performed in the whole (x, y) plane, while the correct integration 
limits are provided by the zero value of the function / outside region ([^). It is instructive 
to specify the direction of by the unite vector ^, which would be shown in a moment to be 
an argument of the Radon transform 



F(^A) = / dt6{t - XxO / d^xe'^'^^fix) =t^xt / dt6{t - xO / d^xe'^"^^ f{x) = 

J —CO J J —CO J 

/oo r _^ roo 

dte'^' d^x5{t-xi)f{x)= dte'^'R{t,0 (12) 
-oo J J~oo 

As a result, the two-dimensional Fourier transform may be represented as a combination of 
one dimensional Fourier transform and Radon transform |^. Consequently, inverting the two- 
dimensional Fourier transform in the polar coordinates and making use of (|12D, the function 
may be expressed through its Radon transform 

1 rco rliT -. roo 

(2^?/. ^"^L *^"'«(P.f) (13) 

This formula already expresses the function f{x,y) in terms of its Radon transform and 
the rest of derivation consists of its simplification. The latter starts with the change of 
variables (p' = (f),p' = p — x^ 

fi^) = 77^.L ^dXj dp'e''^ I d(P'R{p' + ^x,0^ (14) 

x). 



(277)2 7q 



1 rco rco _ 

— XdX dpe'^PR{p, 

271 Jo J-oo 



Here R{p, x) is, obviously, the average of the Radon transform over all the straight lines, 
tangent to the circle of radius p and center in the actual point x. 



*This is actually the core of the relation between DD and SPD, being, respectively, the two-dimensional 
and one-dimensional Fourier transforms of the same matrix element. 
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To present the inversion formula in its final form, one should use the integration by parts 
and substitute the formula for (one-dimensional) Fourier transform of the function sign{X) 

/(f) = _ / sign(\)XdXdpe'^PR(p,x) = — / sign(X)d\dpe'^PR'(p,x) = 

4:71 J-OO ATX J^oo ^ 

r ^R',ip,x) = -- r^(R(p,x)-R{0,x)), (16) 
ZTT J-oo p n Jo p'^ 

where the last equality explores once more the integration by parts and the fact, that 
R{p,x) = R{—p,x). Recalhng the definitions (p|,p!5D, one may express the double distri- 
butions directly in terms of skewed ones (c.f. [|T^): 



dp 

27r^ Jo p2 Jo 



fi^^y) = -7rJlL —L d(j)\cos(l)\{H{p/cos(l) + X + ytgcj), tgcf)) - H{x + ytgcj), tgcj))) 



- / daH{z + x + y^,0-H{x + y^,0)-{l7) 



1 dz 

As one clearly see, to recover the double distribution one indeed should know the skewed 
distribution in an 'extended' region — oo < z,^ < oo. The values ^ < and , \z\, |^| > 1 are 
described by the crossing and symmetry relations. (^,0)- it is easy to check by numerical 
integrations, that GPD obtained from various model DD [jll| may be recovered with the 
reasonable accuracy. 

As it often happened before in the history of the various Radon transform applications 
[^, some of its ingredients were rediscovered in the framework of GPD also. In fact, the 
basic inversion formula in the form of triple integral (similar to ([T3|)) was introduced in 
the framework of GPD by A.V. Radyushkin [l^ (see also ||18|). Also, the polynomiality 



condition |Tl| is well known in the general framework of Radon transform as the Cavalieri 



conditions |I5[. The term is related to the interesting geometric interpretation, so far not 
mentioned in the context of DD. Consider the first [n = 0) moment. It is nothing else 
than the volume of the figure, limited by the surface z = f{x) and the plane x,y. At the 
same time, R{^,p) may be interpreted as a surface of the section of this figure by the plane, 
containing the line p = ^x and z axis. The first moment is then expressing the volume of the 
figure in terms of the surfaces of its intersection with the set of parallel planes. Obviously, 
the direction of these planes is unimportant, so that the moment does not depend on ^. The 
contact with the classical Cavalieri principle may be achieved by the particular choice of 
f{x) = 1 everywhere in the region where it is defined. In that case one is dealing with the 
surface of this region, represented by the lengths of the set of parallel lines. Again it does 
not depend on their direction and ^. 



Polyakov- Weiss terms and two-dimensional magneto- statics 

Let us note that separation of the functions F and G is, generally speaking, impossible, as 
soon as only the information about leading twist SPD (GDA) is available. Indeed, using the 
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integration by parts and assuming that DD turn to zero at the boundary, one can absorb 
the factors z ■ P, z ■ A to the definitions of DD and rewrite in the following way: 

(P'l'*(-|)^-7V{|)|P> = -t j\x 1^^^^^' dye-"'-'-'>'''"N{x,y) (18) 
A,(,,,) = 5^|lJ^)+5^, (19) 

One may absorb the factor 2; ■ P to the definition of the SPD in the analogous way. As a 
result, the 'effective' DD N{x,y) is related by the Radon transform to the 'effective' SPD, 
which is just the derivative ^^qI'^^ of standard SPD. Consequently, it is only possible to 
recover the effective DD, while separation of the individual contributions from F and G is, 
generally speaking, impossible. 

It is interesting, that this ambiguity is completely analogous to the one arising in the 
problem of the recovery of the electromagnetic vector potential from the magnetic field (with 
a zero total fiux). Indeed, denoting G{x, y) = A^, F{x, y) = —Ay, N{x, y) = Bz, we can write 
( p!9D as i? = rotA. Consequently, the effective distribution N{x,y) does not change when F 
and G undergo the following 'gauge transformation' 



G{x,y) G{x,y) - 



dy 
da{x,y) 



dx ' 

a{x,y) = -a{x,-y), (20) 
where the last equation follows from T- (or C-) invariance P). One may try to choose 

a{x,y)= r dtG{t,y), (21) 



so that the G term is completely eliminated. In order to preserve the symmetry properties 
(^ after the gauge transformation, the boundary conditions in ( plf ) should be chosen in a 
following way: 

a{x, y) = \ dtG{t, y) - £ dtG{t, y)^ . (22) 

However, these boundary conditions are the zero ones in the points x = ±(1 — \y\), which in 
this case gurantees the preservation of zero boundary conditions for F and G, only if 

dxG{x,y) = 0, (23) 

s/l-i 

which is indeed true for the case of J = 1 due to the symmetry properties (||). Consequently. 



in this case, when function G is not required by the polynomiality |]TU[, it may be completely 
eliminated by the gauge transformation. At the same time, for / = it is, generally speaking, 
impossible to assume simultaneously the zero boundary conditions at x = 1 — \y\ and x = 
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\y\ — 1 and eliminate G completely. However, it is possible to reduce G{x,y) to the function 
of one variable: 

G{x, y) d{y) = ^ dxGix, y), (24) 

2(1 - \y\) J\yhi 

by the following choice of the gauge: 

a{x, y) = ^ dtG{t, y) - £ dtG{t, y)^ - xd{y) . (25) 

Here the x-independent contribution to G{x, y) is provided by the last term which guarantees 
the zero boundary condition, One can change in this boundary-correcting term the linear 
function to the step one. 



apw{x, y) = - ^ ^ ^ dtG{t, y)~ j dtG{t, y) - sign{x)D{y) 

= e{x<0)r dtG{t,y)-e{x>0) ^'^^ dtG{t,y), (26) 

J\y\-1 Jx 

where D{y) determines the resulting expression for G which now resides at a: = 0: 

G{x, y) 5{x)D{y) = 5{x) dtG{t, y). (27) 



The existence of such a gauge transformation justifies the original suggestion to consider 
this special form of G{x,y). 

Discussion and Conclusions 

We studied crossing properties of Generalized Parton Distributions and showed, that Ra- 
dyushkin's Double Distributions express them in the most simple way. Namely, they simul- 
taneously describe, up to t and W"^ dependence, the DVCS and two-photon production of 
hadron pairs. This allows, in principle, to recover them from the data in the both channels, 
making use of the inverse Radon transform, solving therefore a problem of Radon tomog- 
raphy in the field of non-perturbative QCD. We also considered the ambiguity due to the 
Polyakov- Weiss terms and showed, that it is completely analogous for the gauge ambiguity 
of the vector potential of the static two-dimensional magnetic field. The observable quan- 
tity is in this sense 'effective' DD N{x,y), which may be decomposed to the the F and G 
structures, adopting the particular gauge condition. 

The development of the method may, first of all, include the more refined versions of 
Radon transform (see e.g. numerous Refs. in [ll6[), which would allow to use the limited range 
of angles of Radon projections (limited range of ^). In that case, one may essentially limit 
the required input, and in particular, avoid the use of both channels. Also, one might include 
the possibilities of the singularities of recovered functions [ITTl. Other natural development 



would be to include the twist 3 contributions [20|. Another interseting opportunity may be 



provided by the process 7*7 3n, the GDA |21| for which depend on the two light-cone 
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fractions. The correspondent 'Triple Distributions' would depend on three variables, so to 
recover them one should consider the tomography in three-deimensional space, which (like 
for any odd-dimensional space) is known |]15|, to have the local form due to the Hyugens 
principle. 

Moreover, the Radon transform might be useful in the more general context of quantum 
field theory. Especially interesting seems to be its application to the general problem of 
crossing invariance, as it allows to consider (instead of functions in the different regions of 
their variables) the different projections of the common function, as it was discussed here in 
the particular case of GPD. 

Another interesting application may arise from the fact, that inverse Radon transform 
naturally recovers information about the interior of some region, which makes a contact to 
the famous holographic principle One could add here, that Radon transform technique 



is known for a years to be one of the main tools in the field of optical holography ||16|| . 
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